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ABSTRACT 

This Letter presents the first distance measurement to the massive, semi-detached, eclipsing binary 
LMC-SCl-105, located in the LH 81 association of the Large Magellanic Cloud (LMC). Previously 
determined parameters of the system are combined with new near-infrared photometry and a new 
temperature analysis to constrain the reddening toward the system, and determine a distance of 
50.6 ± 1.6 kpc (corresponding to a distance modulus of 18.52 ± 0.07 mag), in agreement with previous 
eclipsing binary measurements. This is the sixth distance measurement to an eclipsing binary in the 
LMC, although the first to an 0-type system. We thus demonstrate the suitability of 0-type eclipsing 
binaries (EBs) as distance indicators. We suggest using bright, early-type EBs to measure distances 
along different sight lines, as an independent way to map the depth of the LMC and resolve the 
controversy about its three-dimensional structure. 

Subject headings: binaries: eclipsing - stars: distances - distance scale - stars: individual (OGLE 
J053448. 26-694236. 4) - stars: fundamental parameters - galaxies: individual 
(LMC) 



I. INTRODUCTION 

As one of the nearest galaxies to the Milky Way, 
the Large Magellanic Cloud (LMC) has naturally been 
an attractive first rung for the Extragalactic Distance 
Scale. The HST Key Project (jFreedman et al.l [200l 
adopted a distance modulus = 18.50 ± 0.10 mag (cor- 
responding to a distance of 50.1±2.4 kpc) to the LMC, 
which has since become the consensus in the community. 
iSchaefeii (|2008D pointed out that overestimation of error 
bars and band-wagon effects are present in the litera- 
ture, with pre-2001 LMC distance meas urements yield- 
ing values between 18.1 and 18.8 mag (see [Benedict et all 
l2002i) ■ and post-2001 values clustering around the Key 
Project value. Given that different systematic errors 
accompany each method, a careful comparison of the 
distances resulting from different methods is necessary 
to characterize them. Furthermore, there is increasing 
evidence for substantial and complex ver tical structure 
in th e disk of the LMC (see review bv Ivan der Marell 
2006) from studies of red clu mp sta r s (lOlsen fc Salvkl 
2002t iSub ramanian & Subram aniamI |2010[) . Cepheid 
variables ( Nikolaev et alj i2004l ) and RR Lyrae stars 
(jPeicha fc Stanek 2009), which demands further explo- 
ration. 

The only direct, geometrical method available for mea- 
suring distances to stars in the LMC is with eclipsing 
binaries (EBs). In particular, the light curve provides 
the fractional radii of the components, the radial ve- 
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locity semi-amplitudes determine the masses and size of 
the orbit, which together with the effective temperature 
determination (e.g. by comparison with synthetic spec- 
tra), yiel d luminosities and therefore distan ces (see re- 
views by lAndersenI [TOQll: iTorres et al.l [20Tol ). The EB 
distance method has so far been applied to four early- 
B type systems (Guinan et alj 119981: iRibas et all 120021: 
Fitzpatrick et al. 2002, 2003) and one G-type giant sys- 
tem (Pietrzvh ski et al.l i2009) in the LMC, with individ- 
ual uncertainties ranging from 1.2 to 2.2 kpc. Four of 
these systems are located within the bar of the LMC and 
their individual distances are consistent with the quoted 
uncertainties, yielding an error- weighted mean value of 
49.4 ± 1.1 kpc. A fifth system, located several degrees 
away in the north-east quadrant of the disk of the LMC, 
gives a 3cr shorter distance of 43.2 ± 1.8 kpc. 

Figure |T] shows the sp atial distribution of all k nown 
EBs from the OGLE II (iWvrzykowski et al.l [20031) and 
MACHO ffierekas et al.l 120071: IFaccioh et al.il2Ciol mi- 
crolensing surveys of the LMC, and the systems with 
measured distances, overlaid onto the Spitzer SAGE im- 
age in the IRAC 3.6 /im band (Meixner et ah . 200 6). A 
magnitude cut {V < 17 mag) and period cut (> 1.5 
days) were both applied to the EB catalogs to reject 
foreground systems and faint systems whose immedi- 
ate follow up is unrealistic or impossible. The detached 
EBs selected by Michalska & Pigulski (2005) among the 
OGLE II systems as being most suitable for distance de- 
termin ation are a l so sh own. Both the H I kinematic 
center ()Kim et al.l I1998D and the dynamical center (or 
center of the bar; Ivan der Marel etall l200l are over- 
plotted, as is the line of nodes (6 = 129.°9 ± 6.°0deg; 
Ivan der Marel et al. ''2002V 

Motivated by the evidence for vertical structure in the 
LMC and the one discrepant EB distance, we proceed 
to compute the distance to LMC-SC1-10^|. LMC-SCl- 

5 Or OGLE J053448.26-694236.4 = MACHO 81.8881.21 = 
LH 81-72. 
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Fig. 1. — Spatial distribution of known EBs from OGLE II 
and MACHO (blue circles) on the Spitzer 3.6/xm image of the 
LMC. EBs with measured distances are labeled. Yellow circles 
mark the most su itable detached EBs for distance determination 
IIMichalska fc Pigulskii .2005i '): red circles mark the OGLE II bi- 
naries we plan to measure distances to nex t. The H I kine- 
matic center (white "x" ) from lKim et al.l 1199 81') and the dynamical 
center or center of the bar (green "x") from van dcr Marcl ct aLl 
1 2002) are labeled; the sol id line corresponds to the line of nodes 
' van der Marel et al.|[2003 ). Coordinates are given for J2000. 
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105 is a massive, semi-detached, short period {P — 4.25 
days) 0-type system, with component masses of Mi = 
30.9 ± 1.0 Mq, M2 = 13.0 ± 0.7 Mq, and radii of 
Ri = 15.1 ± 0.2 Rq, R2 = 11.9 ± 0.2 Ro (determined 
bv lBonanos|[2009l) . The very accurate measurement of 
the radii (< 2%) renders the system suitable for a dis- 
tance determination, given that EB distances are inde- 
pendent of the usual distance ladder and therefore im- 
portant checks for other methods. However, accurate 
radii are not sufficient for an accurate distance. Accu- 
rate fluxes (i.e. effective temperatures) and extinction es- 
timates are also needed, therefore this Letter sets out 
to determine these quantities and obtain the distance. 
Specifically, Section 2 presents new near-infrared pho- 
tometry of LMC-SCl-105, Section 3 an analysis of the 
spectra with state-of-the-art model atmospheres. Section 
4 the distance determination, and finally. Section 5 a dis- 
cussion of our results. 

2. NEAR-INFRARED DATA 

This study makes use of JHKg observations of 
LMC-SCl-105 obta ined with the CPAPIR camera 
(jArtigau et al.l 120041 ) at the CTIO 1.5-m, as part of a 
synoptic survey of Cepheid variables in the LMC (L. M. 
Macri et al. 2011, in prep.). The EB was observed at 
20 different epochs on 11 nights between 2006 Novem- 
ber 5 and 2007 December 2. Time-series PSF photome- 
try was carried out u sing DAOPHOT and ALLFRAME 
(|Stetsonlll987L 119941 ). Photometric zeropoints were de- 
termined using ~ 2, 500 stars from the 2MASS Point 
Source Catalog, located within 15' of the system and with 
10.5 < Ks < 13.5 mag, while color terms were derived 
using nearly 5 x 10^ 2MASS stars across the entire bar 



Fig. 2.— Phased CPAPIR JHi^^-band light curves of LMC- 
SCl-105. 

of the LMC. Figure [2] shows the calibrated, phased light 
curves of LMC-SCl-105. We calculated error-weighted 
out-of-eclipse mean magnitudes of J = 13.22 ± 0.04, 
H = 13.27 ± 0.04 and Ks = 13.26 ± 0.04 mag. 

3. EFFECTIVE TEMPERATURE ANALYSIS 

An accurate distance measurement to LMC-SC-105 
requires an accurate flux determination for its binary 
components. We proceed to refine the effective tem- 
peratures estimated by Bonanos (2009 j3 with the state- 
of-the-art, NLTE stellar a.tmosphere code F ASTWIND 
(|Santolava-Rev et al.llT997t IPuls et al.l [20051 ). which in- 
cludes the eff'ects of stellar winds and spherical atmo- 
spheric extension. 

The analysis involves a direct comparison between the 
helium lines, which are the main temperature diagnostics 
at these spectral types, plus Ha, to constrain the stellar 
wind, with a complete FASTWIND model grid designed 
to study 0-type stars at the metallicity of the LMC. 
The gri d was developed w ithin the FLAMES-II collabo- 
ration ("Evans et al.ll2010f ) and constructed at the Insti- 
tuto dc Astrofisica de Canarias. Specifically, we derived 
the set of models that provide the lowest x^, using Ha 
and the 10 strongest He I and He II lines availabl^U. The 
synthetic models were downgraded to the instrumental 
resolution of the observed spectra and the projected ro- 

e Tcffi = 35 ± 2.5kK, Toff2 = 32.5 ili 2.5kK, f or log(g) = 3.50 
(fixed), from best fit TLUSTY models I ILanz fc H ubcnv 2003). 

He I AA4026, 4143, 4471, 4713, 4922, 5015, 5875 and He II 
AA4200, 4541, 5411. 
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tational velocities vsini were refined to 160 km s""'^ and 
120 km s~^, for tlie primary and secondary component, 
respectively. We fixed th e surfa ce gravities to the values 
determined by iBonanosI ([20091 ): log(5i) = 3.57 ± 0.02 
and log((72) = 3.40 ± O.Ofl In practice, we rounded the 
values to the first decimal point, to match the 0.1 dex 
step size of the grid. The method provides the stellar 
parameters and their corresponding errors. 

The technique was applied to the two hi ghest S /N 
spectra of LMC-SCl-105 (see lBonanoil2009t ). obtained 
at phases 0.27 and 0.75, i.e. at the first and second 
quadratures. Both phases yielded the same tempera- 
ture for each component, within the errors. Specifi- 
cally, at first quadrature, we found best fit values of 
Teffi = 36100 ± 1000 K, Teff2 = 33200 ± 800 K, while 
at the second quadrature Tosi = 35700 ± 1100 K, 
Toff2 = 33100 ± 900 K. Figures [3] and H show the best 
fit FASTWIND models, plus the effects of the tempera- 
ture errors in the profiles. The synthetic models, which 
only include transitions of H I, He I and He II, provide 
a good match to the observed spectra. Despite not in- 
cluding the Balmer lines in the analysis (except Ha), the 
wings of these lines are in good agreement with the mod- 
els, confirming the accuracy of the \og{g) determination 
from the EB analysis. 
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Fig. 3.— Best fit FASTWIND model (red) of LMC-SCl-105, at 
the first quadrature. The blue (green) lines correspond to models 
with the best fit Teg plus (minus) the la error. The set of lines 
with smaller Doppler shifts corresponds to the primary. 



IBonanosI (|2009[ ) reported changes of the spec tral types 
with phase due to the Struve-Sahade effect (jSticklandi 
Il997f) . the largest being from 07V to 08V for the pri- 
mary, whic h would have an imp act on the temperature of 
~ 2000 K (jMartins et all 120051 ). Our analysis, however, 
does not yield any remarkable differences in temperature 
between the two quadratu res. The reason for this is that 
the classification criteria (jWalborn fc Fitzpatrickiri990[ ) 
hinge on the lines He II A4541, He I A4471, He II A4200, 
and He I-fll A4026, while the FASTWIND analysis av- 
eraged over 10 He I and He II lines in the spectrum. 
The imperfect fits of He II A4200 and He I A4471 by the 
models (see Figure |31), are consistent with a spectral type 
change. 



® Note, the log{g) error bars given in Table 5 of lBonanosI II2009I ') 
incorrectly correspond to the errors in g. 
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Fig. 4. — Same as Figure|31 but for the second quadrature. The 
set of lines with larger Doppler shifts corresponds to the primary. 



At phase 0.75, the secondary star shows important de- 
viations in the cores of H/3 and H7 from the model (see 
Figure 131), which might be due to excess emission arising 
from the slow mass transfer o r the distorted line p rofiles 
of Roche lobe-filling stars fsee lBitner fc Robinson|[2006J . 
Nonetheless, the rest of the He I and He II lines are well 
modeled within the errors. Some of the He II lines (e.g. 
He II A4541) might indicate a slightly higher tempera- 
ture, however these differences lie within the errors. 

4. DISTANCE 

The fiux fx measured at Earth at a certain wavelength 
A from a binary at distance d is given by 



X 10 



-0.4 A{X) 



(1) 



where i?i and R2 are the radii of the two stars and Fi x 
and F2,x the surface fiuxes. The total extinction A (A) 
is a function of the reddening E{B — V), the normalized 
extinction curve k{X - V) = - V)/E{B - V) and 
the ratio of total to selective extinction in the V band, 
Rv = AiV)/E{B -V): 



A{\)^E{B- V) [k {\~V) + Rv] 



(2) 



Having measured the temperatures of the stars from 
the spectra, we computed fluxes and fit to the observed 
magnitudes, using Equation |T] and the best-fit FAST- 
WIND model atmospheres for each quadrature deter- 
mined above. Note that we used the mean radiS of the 
stars instead of their volume radii as better approxima- 
tions to compute their projected surface areaS; 

Foll owing the procedure outlined in Bonanos et al.l 
(|2006D for the detached EB in M33, we calculated syn- 
thetic photometry of the composite spectrum over the 
appropriate Johnso n -Cous ins optical filter fun ctions as 
define d by iBesselll ()1990( ) and calibrated by iLandoltl 
()1992[) . and the 2MASS fiher set. Monochromatic 
fluxes were measured at the isophotal wavelengths (see 
iTokunaga fc Vac"call2005[ ). which best repr esent the flux 
in a passband. We used zeropoints from iBessell et al.l 

^ (rpolo -I- rsido + rback)/3 
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(|1998l Appendix A) and iCohen et all (|2003[ ) to con- 
vert the fluxes to magnitudes. We reddened the model 
spectrum usin g the reddening law parameteriza tion of 
iCardelli et al.l (|1989[ ). as prescribed in S chic ge l et alJ 
([1998'), and simultaneously fit the optica^ and near- 
infrared BVIJHKs photometry. Specifically, we com- 
puted the intrinsic {B — V)q = —0.27 mag from the model 
atmospheres at the isophotal wavelengths, thus yielding 
E{B -V)= 0.11 ± 0.01 mag. 

The value of Rv was determined as the value that min- 
imized the error in the SED fit over the six photometric 
bands. For phase 0.27, we found Rv = 5.8 ± 0.4 and 
for phase 0.75, Rv = 5.7 ± 0.4. The resulting distance 
to LMC-SCl-105 and thus the LMC bar is 50.6 ± 1.6 
kpc (/^ = 18.52 ± 0.07 mag) for the first quadrature and 
50.4 ± 1.6 kpc {p. = 18.51 ± 0.07 mag) for the second 
quadrature. The distances are identical within errors. 
Given the better fit of the FASTWIND models to the 
spectra at first quadrature, we adopt the distance de- 
rived for first quadrature. The fit of the reddened model 
spectrum to the photometry and the residuals of the fit 
are shown in the upper and lower panels of Figure [5l 
respectively. The error in the distances was computed 
by a bootstrap resampling procedure. We repeated the 
spectral energy distribution (SED) fitting procedure 1000 
times for each quadrature, by randomly selecting (using 
Gaussian sampling) all the parameters within their er- 
rors. We adopt the a of the resulting Gaussian distribu- 
tion as the uncertainty in the distance. 

We tested the robustness of our reddening and distance 
results, by first fitting the BVI photometry alone, which 
yielded an identical value for the distance (50.8 ±1.6 
kpc or /X = 18.53 ± 0.07 mag, with Rv = 5.7 ± 0.4), 
thus demonstrating the consistency of the near-infrared 
with the optical photometry. Next, if we fix Rv = 3.1, 
the best fit value for E{B — V) — 0.18 mag, resulting 
in a distance of 51.9 ± 1.6 kpc (/x = 18.58 ± 0.07 mag), 
i.e. in agreement with our reported result, within errors. 
If instead we assume Rv = 3.1 and fix E{B — V) = 
0.11 ±0.01 mag (based on our photometry and the model 
spectraj3, we would derive a much larger distance of 
55.2 kpc (^ = 18.71 mag), which yields a SED fit error 
of 0.05 mag (versus 0.01 mag) that is inconsistent with 
the photometry. The validity and implications of the 
high value of Rv that we have measured are discussed in 
the following Section. 

The error quoted above for R y was estimated usin g 
the Bayesian code CHORIZOS (jMai'z-Apellanid 120041 ). 
The available BVIJHKs photometry was given as input, 
with Teff in the range 33000-36000 K and log(g) fixed 
to 3.50, from TLUSTY models. The code yielded best fit 
mean values (for a single star) of Tog = 34500 ± 1100 K, 
-RA5495 = 5.4 ± 0.4 and £;(A4405 - A5495) = 0.10 ± 0.01 
mag, consistent with the values we derived. 

5. DISCUSSION 

LMC-SCl-105 is located in the LH 81 association 
(|Massev et al.ll2000D . near the center of the LMC bar. 
It contains two early 0-type stars and three Wolf-Rayet 

Bmax = 12.81 ± 0.01 mag, Vmax = 12.97 ± 0.01 mag, /^ax = 
13.04 ±0.01 mag tWvrzvkowski ct al. 2003). 

Note, our E{B — V) value is consistent with the range (0.13- 
0.23 mag) measured bv lMassev et al.l l|2000l ) for 34 stars in LH 81. 
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Fig. 5. — Upper panel: fit of the reddened EB model spectrum 
(for phase 0.27) to the BVIJHKs photometry. Lower panel: resid- 
uals of the SED fit, in terms of the flux ratio. Error bars correspond 
to the photometric error for each band in flux units. The best fit 
values of E{B - V) = 0.11 ± 0.01 mag and Ry = 5.8 ± 0.4 yield a 
distance modulus to the EB and thus the LMC bar of 50.6 it 1.6 
kpc {^l = 18.52 ± 0.07 mag). 

systems, one of which was recently found to be an EB 
(jSzczvgiel et al.l[2"010l ) . Furthe rmore, this as sociation re- 
sides in the superbu bble N 154 ()Henizelll956D = DEM 246 
(|Davies et al.l 119761 ). We have determined a large value 
of Rv = 5.8 ± 0.4 toward LMC-S Cl-105, howeve r such 
high values are not uncommon. ICardelh et aP (fl989l ) 
find 5 < Rv < 5.6 for 6 out of the 29 OB stars in their 
sample, while lFitzpatrick fc Massal (f2007) find i?^/ > 5 
for 12 out of the 328 stars in their sample. Large val- 
ues of Rv simply imply larger dust grain sizes, which 
are expected to occur in dense regions of the interstellar 
medium due to accretion and coagulation of grains. We 
therefore conclude that the environment in which LMC- 
SCl-105 resides has large dust grains. 

In this Letter, we have determined the distance to 
LMC-SCl-105 and consequently the LMC bar to be 
50.6 ± 1.6 kpc (/Lt = 18.52 ± 0.07 mag). The agreement 
we find with previous EB distances to systems in the bar 
with different spectral types testifies to the robustness 
of the EB method and its potential as a powerful, in- 
dependent distance indicator. Furthermore, it confirms 
that 0-type (and semi-detached) EBs are suitable for 
distance determination, i.e. that the fluxes predicted by 
FASTWIND are indeed accurate. EB-based distance de- 
term inations to M31 (iRibas et al.1 120051 : iVilardell et~all 
HMQ) and M33 (jBonanos et al.l l2006') can therefore pro- 
vide an independent absolute calibration of the Extra- 
galactic Distance Scale. Future distance determinations 
to EBs in the LMC (e.g. those marked in Figure [Ij, will 
additionally provide Rv values in different environments 
of the LMC. Finally, we suggest using bright, early- type 
EBs to measure distances along different sight lines to 
the LMC, as an independent way to map its depth and 
resolve the controversy about its vertical structure. 
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